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Is Fermi-surface nesting the origin of superconductivity in iron pnictides?: 
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We study whether Fermi-surface (FS) nesting can give rise to high-temperature superconduc- 
tivity in iron pnictides. Starting with ab initio construction of an effective four-orbital model, we 
employ the fluctuation-exchange approximation to show that FS does not necessarily favor the 
stripe antiferromagnetic order observed in experiments, especially for realistic electronic cor- 
relations. If superconductivity in iron pnictides is magnetically mediated and has fully-gapped 
sign-reversing s-wave symmetry, our results suggest that the pairing interaction does not arise 
only from FS nesting and exchange interactions between local moments in the Fe 3d orbitals 
may play a crucial role. 
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The seminal discovery of superconductivity in 
LaFeAsO doped with fluorine^ and subsequent updat- 
ings of the record of the transition temperature (Tc) in 
related iron pnictides and iron chalcogenites^ have at- 
tracted great interest in these compounds. Since there is 
close proximity of antiferromagnetic (AF) order and su- 
perconductivity in the phase diagram, the possibihty of 
magnetically mediated superconductivity has been inten- 
sively studied.'^ 

When we consider the interplay between magnetism 
and superconductivity, one of the key issues to be clar- 
ified is, whether the magnetism should be understood 
in terms of locahzed spins or itinerant electrons. While 
there have been a variety of proposals and calculations 
for the so-called J1-J2 model based on the the strong cou- 
pling scenario,"* it has also been claimed that the weak 
coupling approach should work better, ^"^* and the prob- 
lem which picture describes the system more properly 
has been an issue of hot debates.'^ 

Among many studies based on the weak coupling sce- 
nario, recently, a detailed calculation by the random 
phase approximation (RPA) for a five-orbital model de- 
rived from ab initio calculation based on the local density 
approximation (LDA) has been performed. There, it has 
been discussed how changes in the crystal structure (i.e., 
the atomic configuration) affects superconductivity,^ and 
it has been concluded that the weak coupling approach 
succeeds to explain the qualitative tendency of pairing 
instability in iron pnictides: when the pnictogen height 
is low (high) as in the case of LaFePO (LaFeAsO), low 
(high) Tc nodal (nodeless) pairing is favored. 

On the other hand, estimates of the interaction pa- 
rameters by various ah initio methods suggest that iron 
pnictides are not weakly but moderately correlated. ^'^""'^^ 
Indeed, photoemission spectrum indicates that the mass 
enhancement due to the electron correlations is as large 
as 2 ~ 3,^^ whose effect is neglected in RPA. Since Fermi- 
surface (FS) nesting plays a crucial role in the weak cou- 



pling approach, and the orbital/momentum dependence 
of the self-energy correction can seriously affect the nest- 
ing condition, it is of great interest to see whether the 
weak coupling approach based on FS nesting survives 
even if we consider the self-energy correction. 

One of the standard ways to investigate the ef- 
fect of the self-energy correction is to employ the dy- 
namical mean field theory (DMFT).!"^ While several 
DMFT calculations combined with LDA have been per- 
formed,^^' the recent calculation with realistic in- 
teraction parameters'^ for LaFeAsO has shown that the 
orbital-dependent DMFT correction lowers the level of 
d^2_y2 at the F point in the Brillouin zone (hereafter, the 
local coordinate system centered at Fe is such that the 
X- and y-axis point towards neighboring As). As is dis- 
cussed below, in such a situation, the nesting condition 
is expected to be worse, which suggests that the nesting 
scenario is not necessarily robust against the self-energy 
correction. 

While the self energy in DMFT is momentum indepen- 
dent, the momentum dependence of the self energy is also 
crucial for FS nesting. The purpose of the present study 
is to investigate the effect of momentum-dependent self 
energy for realistic band structure and electronic corre- 
lations by means of the fluctuation-exchange approxima- 
tion(FLEX).'^ Although FLEX can become problematic 
to address some high energy features (e.g., it fails to de- 
scribe a Mott insulator), FLEX is generally expected to 
be reliable for description of low energy physics origi- 
nating from FS nesting, so that this tool is expected 
to work for the examination whether superconductivity 
in iron pnictides, which are moderately correlated, can 
indeed be understood in terms only of FS nesting. 

Let us move on to the construction of the effective 
model which we employ in the following calculation. As 
for FLEX calculation for iron pnictides, while there have 
already been several studies for simplified models,^ one 
of the present authors recently performed a calculation 
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for a five-orbital model based on LDA.* There, he set 
the intra-orbital Coulomb interaction U and the Hund 
coupling J to be 1.44 eV and 0.24 eV, which are rela- 
tively small compared to those estimated by an ab initio 
method. ^^'^^ In fact, it has been recognized that FLEX 
becomes problematic for the five-orbital model in the 
moderately correlated (realistic) regime since some un- 
controllable treatment is needed to get a reasonable spec- 
tral function. Namely, we have to introduce an artificial 
level shift for 1^22 to the original one-body Hamiltonian,^ 
otherwise the dz^ orbital makes a large FS around the F 
point which is not consistent with the experiments.^^ 

The reason why we have to shift the level of can 
be understood as follows. Since the the d^^ orbital in the 
five-orbital model is localized and correlated^^ (due to 
weak hybridization with As Ap), the correlation already 
taken in LDA is expected to be large for ^^2 . Thus, while 
we have to introduce the so-called double counting (DC) 
term^^ to combine LDA with FLEX, the DC term of ^^2 
should be larger than those of dy^ , d^z and d^2_y2 (which 
make the FS in LDA). However, it is a non-trivial task 
to evaluate the orbital dependent DC term from first 
principles for the five-orbital model, especially when the 
system is moderately correlated. 

In order to avoid this subtle problem and make our 
analysis on the validity of the nesting scenario be clearer, 
in the present study, we first construct an effective four- 
orbital model for d^y, dx2.y2, dyz, and dxz- As we can 
see in Fig. 1, since dz^ basically lies below the Fermi 
level, it is expected to be irrelevant for the physics of the 
FS nesting. In actual model construction, as was done in 
Ref.,^ we obtain the band structure of LaFeAsO with the 
Quantum-ESPRESSO package, and made maximally 
localized Wannier functions^^ for the four orbitals to rep- 
resent the Kohn-Sham Hamiltonian in terms of these or- 
bitals. 

In table I, we list the transfer integrals in the one- 
body part of the Hamiltonian. The onsite energies of 
dyzidxz), dx2.y2, dxy s^YC 0.09, 0.28, and O.lleV with re- 
spect to the Fermi level. The spread of Wannier func- 
tions (V(r2) - (r)2) are 2.16, 2.55 and 2.11 A, while in 
the five-orbital model, they are 2.05, 2.36, and 1.73 A, 
respectively. 

In Fig.l, we plot the band dispersion (in the unfolded 
Brillouin zone^ for the unit cell which contains only one 
Fe atom) of the four-orbital model, comparing with that 
of the five-orbital model. We can see that the FS is 
almost the same as that of the original five-orbital model. 
Thus we may expect that if the nesting scenario indeed 
captures the essential feature of the pairing mechanism, 
the four-orbital model should also be able to describe 
superconductivity in iron pnictides. 

For the many-body part of the Hamiltonian, we con- 
sider the standard interactions that comprise the inter- 
and intra-orbital Coulomb interaction {U, U') and the 
Hund coupling and the pair-hopping (J, J'). We assume 
that U = U'+2J and J = J' . Following the standard pro- 
cedures of multi-orbital FLEX,^ we solve the linearized 
Eliashberg equation, and obtain the gap function along 
with the eigenvalue A. Tc corresponds to the temperature 
where A reaches unity. 32x32x4 three-dimensional k- 
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Fig. 1. (a)Band dispersion of the four-orbital model (solid line) 
and the five-orbital model (dotted line) for LaFeAsO. (b)Fermi 
surface of the five-orbital model labeled as ai, 02, Pi, f32 and 
7, together with nesting vectors qi and q2. (c)A plot similar to 
(b), where the d^2_j^2 -character is dominant is shown by dotted 
lines. 



point meshes and 1024 Matsubara frequencies are taken, 
while the temperature (T) is set to be 0.02 eV. Hereafter 
we denote the spin susceptibility, the gap function, and 
the squared absolute value of the Green function for the 
lowest Matsubara frequency (ia;„ = for the spin sus- 
ceptibility and inT for the latter two) and fc^ = as 
Xsikx,ky), A{kx,ky), and g{kx,ky), respectively. 

First, we compare the four-orbital model and the five- 
orbital model by RPA. Both models have a peak in the 
spin susceptibility which corresponds to stripe AF insta- 
bility, and the fully-gapped s-wave with sign reversing 
(the so-called s-^-wave) is most dominant for the un- 
doped case (electron-filling n is 4.0 for the four-orbital 
model and 6.0 for the five-orbital model). Interaction pa- 
rameters are taken as, U — 1.2, U' — 0.9, and J = 0.15 
eV. In Fig. 2, we plot A of the the four(five)-orbital 
model, for the 1st (2nd) to 3rd (4th) bands from bot- 
tom to top in the band representation. We see that the 
gap functions of these two models are similar to each 
other in RPA. 

Next let us move on to FLEX. In Fig. 3, we com- 
pare RPA and FLEX results of g (which takes large val- 
ues at FS, see also Fig. 1(c)) for dx2_y2 in the orbital 
representation and Xs- Here we set U = 1.2, U' = 0.9 
and J = 0.15 eV. Although we have a sharp peak in 
Xs at (kx,ky) = (tt, 0) and (0, tt) in RPA, these peaks 
are severely suppressed in FLEX and a broad structure 
around (tt, tt) becomes dominant. While the fully-gapped 
s±-wave is dominant in both RPA and FLEX, A is much 
smaller in FLEX(=0.18) than in RPA(=0.83). These re- 
sults suggest that, in order to have strong magnetic flue- 
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Table I. Hopping integrals t(AX, Ay; ^, v) in units of eV. [AX, AY] denotes the in-plain hopping vector (note that the X- and K-axis 
point towards neighboring Fe atoms), and (/i, the orbitals. Uy, /, and (T^ corresponds to t(AX, — AY; ^, i^), t(— AX, — AY; /i, !^), 
t(AY, AX; /i, f), respectively, where '±' and '±(/i',v')' in the row of {fi,!^) mean that the corresponding hopping is equal to 
±t{AX, AY; fi, y) and ibt(AX, AY; u'), respectively. This table, combined with the relation t(AX, AY; ^, v) = t{-AX, -AY; u, /i), 
gives all the in-plain hoppings > 0.02eV up to the fifth neighbors. 
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Fig. 2. Gap functions for the four-orbital model (a) (the five- 
orbital model (b)) of LaFeAsO, for the 1st (2nd) to 3rd (4th) 
bands from bottom to top in the band representation with n = 
4.0 (6.0). Solid lines represent the Fermi surface. 

tuations and pairing instability from the nesting between 
small Fermi pockets, generally we need sufficiently large 
interactions. 

In Fig. 4, we show the result of FLEX for U — 1.8, 
U' = 1.2 and J = 0.3. In contrast with Fig. 3(b), we have 
a sharp peak in Xs (Fig. 4(a)). However, the position of 
the peak is located at (tt, tt) which corresponds to the G- 
type AF instability. Here, rather than the fully-gapped 
s_l--wave pairing, the so-called gapless s±-wave pairing 
(whose |A| becomes small on the Fermi surface^) and a 
d-wave pairing are more favored (A is ^ 0.34 for these 
pairings, and see Fig. 4(c) in which A of these pairings 
for the third band in the band representation are plot- 
ted). For larger interaction parameters, the G-type AF 
instability is more and more enhanced and the gapful 
s±-wave pairing never becomes dominant. 

These results can be understood as follows. For un- 
doped LaFeAsO, as we see in Fig. 1(b) and (c), there are 
four kinds of Fermi pockets: two hole pockets, ai and 
a2, around (0,0) and electron pockets, /3i, /32, around 




Fig. 3. RPA(a) and FLEX(b) results for the squared absolute 
value of the Green function for the d^2_y2 orbital in the orbital 
representation (left) and spin susceptibility at the lowest Mat- 
subara frequency (right) for the four-orbital model with U = 1.2, 
U' = 0.9, and J = 0.15 eV. 

(tt, 0) and (0, tt) as well as a hole pocket 7 around (tt, tt) 
(which disappears by ~ 10% doping of electrons). While 
a is made from c?j,zand d^z, P is made from dyz, dxz and 
c?^2_j,2, and 7 is made from d.j.2_y2 (Fig. 1(c)). On top 
of the scattering channels which connect a and /3 or 7 
and f3 (qi and q2 in Fig. 1(b)), there is also a channel 
between dx2_y2 in /3i and (32 (qa in Fig. 1(c)). As is dis- 
cussed in Ref.,^ while q2 is the key to induce the strong 
stripe AF fluctuations, in the absence of 7, qs suppresses 
the stripe AF fluctuations and the fully-gapped s-|--wave 
pairing. In fact, the low energy electronic structure in 
FLEX is quite similar to that in RPA: the spectral func- 
tions obtained by these approximations have similar size 
of structures in the Brillouin zone around (7r,0), (0,7r), 
(0, 0) and (tt, tt) for low frequencies. However, the self- 
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energy effect for d,j.2_y2 around (tt, tt) (where tlie LDA 
band dispersion is flat) is significantly large and severely 
suppresses q2, which is dominant in RPA. While effects 
which have not been considered (such as spin-lattice cou- 
pling, off-site Coulomb interactions, etc.) may change the 
situation, the present study suggests that the nesting sce- 
nario is not necessarily robust against the momentum- 
dependent self-energy correction. 




Fig. 4. FLEX results for the four-orbital model with U = 1.8, 
U' = 1.2, and J = 0.3 eV. Spin susceptibility (a), squared ab- 
solute value of the Green function for the d^2_y2 orbital (b) and 
the gap function of the d^2_y2 orbital for the s±- and d-wave 
pairing at the lowest Matsubara frequency (c) are shown. 

To summarize, we study whether FS nesting can be the 
origin of spin-fluctuation-mediated superconductivity in 
iron-pnictides by means of FLEX. Our result suggests 
that the nesting scenario has a serious dilemma: since 
all the FS pockets are small, we need sufficiently strong 
interactions to have large magnetic fluctuations from the 
FS nesting. However, due to the self-energy corrections, 
the contribution of the low-energy d,j.2_y2 state at (tt, tt) 
which is the key to make dominant stripe AF fluctuations 
in the weak coupling regime is severely suppressed for 
the moderately correlated regime. The remaining dx2.y2 
states at (0, tt) and (tt, 0) induce the G-type AF fluctua- 
tions which favor the d- or nodal s-wave pairing and the 
high Tc nodeless s±-pairing is difficult to arise. 

Acknowledgment 

We would like to thank H. Aoki, P. Hansmann, M. 
Imada, G. Sangiovanni, A. Toschi for fruitful discussions 
and/or critical reading of the manuscript. This work is 
supported by a Grant-in- Aid for Scientific Research on 
Priority Areas "Super Clean Materials" under Grant No. 
20029014 from MEXT, Japan. Numerical calculations 
were performed at the facilities of the Supercomputer 
center, ISSP, University of Tokyo. 



1) Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130 (2008) 3296 

2) for review, see e.g., K. Ishida, Y. Nakai and H. Hosono, J. 
Phys. Soc. Jpn 78 (2008) 062001. 

3) L L Mazin and J. Schmalian, Physica C 469(2009) 614. 

4) Q. Si and E. Abrahams, Phys. Rev. Lett. 101 (2008) 076401, 

C. Fang, H. Yao, W-F. Tsai, J. Hu, S. A. Kivelson, Phys. 
Rev. B 77 (2008) 224509, C. Xu, M. MuUer and S. Sachdev, 
78 (2008) 020501(R), K. Seo, B. A. Bernevig, and J. Hu, Phys. 
Rev. Lett. 101 (2008) 206404. 

5) K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. Kontani, 
and H. Aoki, Phys. Rev. Lett. 101 (2008) 087004, erratum: 
Phys. Rev. Lett. 102 (2009) 109902(E). 

6) K. Kuroki, H. Usui, S. Onari, R. Arita, and H. Aoki, Phys. 
Rev. B 79 (2009) 224511. 

7) S. Graser, T. A. Maier, P. J. Hirschfcld, and D. J. Scalapino, 
New J. Phys. 11 (2009) 025016. 

8) H. Ikeda, J. Phys. Soc. Jpn. 77 (2008) 123707, erratum: J. 
Phys. Soc. Jpn. 78 (2009) 038001 

9) Z. J. Yao, J. X. Li and Z. D. Wang, New. J. Phys. 11 (2009) 
025009, R. Sknepnek, Phys. Rev. B, 79 (2009) 054511. 

10) A. V. Chubukov, D. Efremov, L Eremin, Physical Review B 
78 (2008) 134512, F. Wang, H. Zhai, Y. Ran, A. Vishwanath, 

D. -H. Lee, Phys. Rev. Lett. 102 (2009) 047005, F. Wang, H. 
Zhai, D.-H. Lee, Europhysics Lett. 85 (2009) 37005. 

11) T. Nomura, J. Phys. Soc. Jpn. 78 (2008) 034716. 

12) for recent experiments which suggest the existence of nodes 
in the pairing gap function of LaFePO, see M. Yamashita, N. 
Nakata, Y. Senshu, S. Tonegawa, K. Ikada, K. Hashimoto, H. 
Sugawara, T. Shibauchi, Y. Matsuda, arXiv: 0906.0622. J. D. 
Fletcher, A. Serafin, L. Malone, J. Analytis, J-H Chu, A. S. 
Erickson, I. R. Fisher, and A. Carrington, Phys. Rev. Lett. 102 
(2009) 147001. C. W. Hicks, T. M. Lippman, M. E. Huber, J. 

G. Analytis, J.-H. Chu, A. S. Erickson, L R. Fisher, and K. 
A. Moler, arXiv: 0903.5260 

13) K. Nakamura, R. Arita, M. Imada, J. Phys. Soc. Jpn. 77 (2008) 
093711. 

14) T. Miyake, L. Pourovskii, V. Vildosola, S. Biermann, and A. 
Georges, J. Phys. Soc. Jpn. 77 (2008) Supplement C 99. 

15) V. I. Anisimov, Dm. M. Korotin, M. A. Korotin, A. V. 
Kozhevnikov, J. Kunes A. O. Shorikov, S. L. Skornyakov, and 
S. V. Streltsov, J. Phys. Cond. Matt 21 (2008) 075602, V. L 
Anisimov, E. Z. Kurmaev, A. Moewes, I. A. Izyumov, Physica 
C 469 (2009) 442. 

16) see e.g., D. H. Lu, M. Yi, S.-K. Mo, J. G. Analytis, J.-H. Chu, 
A. S. Erickson, D. J. Singh, Z. Hussain, T. H. Geballe, I. R. 
Fisher, Z.-X. Shen, Physica C 469 (2009) 452, W. Malaeb, 
T. Yoshida, T. Kataoka, A. Fujimori, M. Kubota, K. Ono, H. 
Usui, K. Kuroki, R. Arita, H. Aoki, Y. Kamihara, M. Hirano, 

H. Hosono, J. Phys. Soc. Jpn. 77 (2008) 093714. 

17) W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62 (1989) 324; 
A. Georges and G. Kotliar, Phys. Rev. B, 45 (1992) 6479; A. 
Georges W. Krauth G. Kotliar and M. Rozenberg, Rev. Mod. 
Phys. 68 (1996) 13. 

18) K. Haule, J. H. Shim, and G. Kotliar, Phys. Rev. Lett. 100 
(2008) 226402. 

19) S. L. Skornyakov, A. V. Efremov, N. A. Skorikov, M. A. Ko- 
rotin, Yu. A. Izyumov, V. I. Anisimov, A. V. Kozhevnikov, 
and D. Vollhardt, arXiv:0906.3218 

20) M. Aichhorn, L. Pourovskii, V. Vildosola, M. Ferrero, O. Par- 
collet, T. Miyake, A. Georges, and S. Biermann, Phys. Rev. 
B 80 (2009) 085101 

21) N. E. Bickers and D. J. Scalapino, Ann. Phys. (N.Y.) 193 
(1989) 206. 

22) It should be noted that FLEX can become problematic even 
in the weak coupling regime. For example, it does not work 
well for materials such as Sr2Ru04 where FS nesting is not 
good, see e.g., Y. Yanase, T. Jujo, T. Nomura, H. Ikeda, T. 
Hotta and K. Yamada, Physics Reports, 387 (2003) 1. 

23) C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. Ratcliff II, J. 
L. Zarestky, H. A. Mook, G. F. Chen, J. L. Luo, N. L. Wang, 
and P. Dai, Nature (London) 453 (2008) 899. 

24) V. I. Anisimov, F. Aryasetiawan and A. I. Lichtenstein, J. 



J. Phys. Soc. Jpn. 



Letter 



Author Name 5 



Phys. Condensed Matter 9 (1997) 767 

25) S. Baroni, A. Dal Corso, S. de Gironcoh and P. Giannozzi, 
http://www.pwscf.org/. Here we adopt the exchange correla- 
tion functional introduced by J. P. Perdew et al.[Phys. Rev. 
B 54 (1996) 16533], and the wave functions are expanded by 
plane waves up to a cutoff energy of 40 Ry. 10^ fe-point meshes 
are used with the special points technique by H. J. Monkhorst 
and J. D. Pack [Phys. Rev. B 13 (1976) 5188]. 

26) N. Marzari and D. Vanderbilt, Phys. Rev. B 56 (1997) 12847; I. 



Souza, N. Marzari and D. Vanderbilt, Phys. Rev. B 65 (2002) 
035109. In this study, we used the code developed by A. A. 
Mostofi, J. R. Yates, Y.-S. Lee, L Souza, D. Vanderbilt and N. 
Marzari (http://www.wannier.org/) to construct the Wannier 
functions. 

27) V. Vildosola L. Pourovskii, R. Arita, S. Biermann, A. Georges 
Phys. Rev. B., 78 (2008) 064518. 

28) It should be noted that the band dispersion of the five-orbital 
model is exactly the same as that of ab-initio calculation.^' 



